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Abstract 

We study the effects of the finite size of baryons on the equa¬ 
tion of state of homogeneous hadronic matter. The finite extension 
of hadrons is introduced in order to improve the performance of field 
theoretical models at very high densities. We simulate the in-medium 
averaged baryon-baryon strong repulsion at very short distances by 
introducing a Van der Waals like normalization of the fields. This is 
done in the framework of the Quark Meson Coupling model that al¬ 
lows to take care of the quark structure of baryons. Since within this 
model the confinement volume evolves with the fields configuration, 
the treatment is not equivalent to a simple hard-core potential. We 
investigate the phase transition to quark matter and the structure of 
neutron stars. We have found significant corrections at high densities. 
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1 Introduction 


Investigation of hadronic matter at extreme conditions of density and tem- 
peratnre is a current issue of research, since the study of its properties will 
eventually shed some light over the recovering of QCD symmetries PIE]- The 
phase diagram of hadronic matter is expected to be very complex, exhibit¬ 
ing exotic phases like superfluidity, meson condensates, dibaryon condensate, 
etc. The gradual emergence of quark droplets would hnally lead to a tran¬ 
sition to deconhned quark matter. Every one of these phenomena affect the 
equation of state (EOS) and could have macroscopic manifestations as, for 
example, in the structure of stars. 

In most hadronic matter studies baryons are assumed to be point-like. 
This can be justihed because at densities below the nuclear matter saturation 
density the hnite volume effects are expected to be small. However, the 
spatial extension of baryons was recognized as an essential point in the study 
of the collective phenomena at very high densities jS]. 

It is worthy to mention that there are only a few field theoretical models 
which consistently include the baryonic spatial extension. The most com¬ 
monly used are the Skyrme and bag-like models. In the first case the inclusion 
of hnite baryon density effects is not straightforward due to the topological 
character of its solutions |3]. On the other hand, further rehnements of the 
original MIT bag model allow to deal with medium effects upon the hadron 
structure 000 Eli. Within this scheme there were recent efforts to include 
the repulsion between overlapping bags, by means of effective short-ranged 
quark-quark correlations 0. 

The authors attempted to take into account hnite volume correction pre¬ 
viously, using a Van der Waals-like method to study nuclear matter with 
lambda hyperons PI- The total volume appearing in the thermodynamical 
quantities was replaced by the available volume, in accordance with related 
investigations nulla ESI EH iini 

In the present work we generalize this approach to study the properties 
of hadronic matter including the octet of low lying baryons, and to check 
out their inhuence over the transition to quark matter. The resulting EOS 
is applied to study the composition of neutron stars. Following our previous 
study mi. we introduce these corrections at the level of the normalization 
of the baryon helds in the Quark Meson Coupling (QMC) model 00. The 
motivation for such a procedure is to parameterize in compact form the 
strong baryon-baryon repulsion at very short distances. We focus on the 
high density regime, thus it is justihed to consider the statistical average of 
the interaction instead of looking at its details. 

In the QMC model the size of the conhning volume has its own dynam¬ 
ical evolution which takes into account the baryonic density and the helds 
conhguration, and it is obtained in a self-consistent calculation. Thus it 
can somehow be regarded as an ehective degree of freedom. The bag radius 
changes smoothly with the medium properties, and in this sense the hnite 
volume corrections introduced here cannot be considered as a hard core in¬ 
teraction. 

In the next section we give a resume of the QMC model and we introduce 


the excluded volume correlations. In section 3 we describe neutron star 
matter and the phase transition to quark matter. Numerical results and 
discussions are given in section 4, and conclusions are drawn in section 5. 


2 The Quark Meson Coupling Model 

Relativistic Hadron Field Theories provide a good description of nuclear mat¬ 
ter near the saturation density, and of hnite nuclei too. For this purpose only 
a small number of free parameters is required. Mean field approximation 
(MFA) is suited for the aim of these theories, since the treatment of matter 
at medium and high densities do not require the detailed structure of the 
interactions. Within this framework the QMC model PIE] can be viewed as 
an extension of the Quantum Hadro-Dynamics (QHD) models [THl ITT] . 

In the QMC model baryons are represented as non-overlapping spherical 
bags containing three valence quarks; the bag radius changes dynamically 
with the medium density. Baryons interact by the exchange of a, lo and p 
mesons coupled directly to the conhned quarks. It has been found that these 
extra degrees of freedom, provided by the internal structure of the baryon, 
lead to quite acceptable values of the nuclear matter compressibility at sat¬ 
uration. Despite of the explicit quark helds in the QMC model, hadronic 
thermodynamical properties are evaluated in such a way that baryons are 
handled as point-like particles with an effective mass which depends on 
the a held. 

In the MFA the Dirac equation for a quark of havor g, (g = u,d,s), of 
current mass rriq and Jg third isospin component, is given by 


(n'-S, - - syi iX - = 0. 


( 2 . 1 ) 


In this equation all meson helds have been replaced by their mean held values. 
Mesons couple linearly only to non-strange quarks, i.e. = Pp = 0. 

Therefore the parameters rriq* are given by 




( 2 . 2 ) 




For a spherically symmetric bag of radius Rh representing a baryon of class 
b, the normalized quark wave function T^(r, t) is given by 



(2.3) 


where is the quark spinor and 



(2.4) 



(2.5) 




(2.6) 







with flqb = [Xqf^ + The eigenvalue Xqb is solution of the equation 


joiXqb) Pqjl{Xqb), ( 2 - 7 ) 

which arises from the boundary condition at the bag surface. 

In this model the ground state bag energy is identihed with the baryon mass 

MS, 


MS = 


Kb o 


( 2 . 8 ) 


where is the number of quarks of flavor q inside the bag. The bag constant 
B represents the difference of energy per unit volume between the vacuum 
with and without broken QCD symmetry. It is numerically adjusted to get 
dehnite values for the proton bag radius. The zero-point motion parameters 
ZQb are hxed to reproduce the baryon spectrum at zero density. 


Eq. (j2.8j) shows that the baryon effective mass is a function of the bag 
radius Rb- In the original MIT bag calculations Rb is a constant hxed at 
zero baryonic density, but in the QMC it is a variable dynamically adjusted 
to reach the equilibrium of the bag in the dense hadronic medium. We use 
the equilibrium condition proposed in ref. HU, which results by imposing 
a vanishing net hux of the energy-momentum tensor through the surface of 
the bag immersed in the dense hadronic medium. This gives 


1 (dMb*\ 

Ml V dR, 



dkk^ 


(2.9) 


where ^ = 1. The interested reader can hnd the detailed derivation of this 
relation in HU. 

This result rehects the balance of the internal pressure of the bag with 
the baryonic contribution to the total external pressure, represented by the 
left and right sides of Eq. dzni), respectively. 

The factor ^ will be redehned below when excluded volume effects will be 
considered. 

It must be noted that Eq. (E3D differs from the standard QMC condition 

[HI, i.e. 



Both prescriptions coincide only in the case of vanishing density. Ea. (12.9|1 could 
be more appropriate for dealing with hnite density calculations, nevertheless 
there remains to elucidate the problem of overlapping bags as the density 
grows. Therefore it becomes necessary to introduce some further considera¬ 
tions into the formalism to handle this feature. 

Once Mb* has been dehned microscopically, the hadronic thermodynamics 
in the QMC model resembles that of the Quantum Hadrodynamics. In the 
MFA for homogeneous inhnite static matter all meson helds are replaced by 
their averaged values, i. e. 










( 2 . 12 ) 


( 2 . 10 ) 


( 2 . 11 ) 


where a = 1, 2, 3 runs over all isospin directions and /| is the third isospin 
component of baryon b. In our calculations we have used the values = 
550MeV^, = 783MeV", and rrip = 770MeV for the meson masses. 

The dispersion relation for the b-baryon is 



for particle (+) and antiparticle (—) solutions. Within the MFA at zero tem¬ 
perature only the particle solutions contribute. 


The scalar (n^) and baryonic (n^) densities are dehned with respect to 
the ground state of the hadronic matter \GS > composed of baryons hlling 
the Fermi sea up to the state with momentum kb 



(2.14) 


(2.15) 


In Eqs. (^rm and (EH) the factor 7} is included for future use and it takes 
the value = 1 for point-like baryons. 

In the next section we describe hadronic matter in /^-equilibrium and 
electrically neutral, therefore we also consider leptons treated as free Dirac 
particles. The leptonic density ni is related to the Fermi momentum ki by 


n' =< >= 


(2.16) 


Given a distribution of baryonic species we can calculate the total energy 
density en and pressure Pq of hadronic matter for point-like baryons 



(2.17) 


a = 1 *’"’ - ch, 


b I 


(2.18) 













where /Iq = k^ikb) (see Eq. ()2.13j) 'l and /i* are the chemical potentials for 
point-like baryons and leptons, respectively. 

The QMC model has been widely used to describe the dense hadronic 
matter, it is based in a great extent on the assumption of non-overlapping 
bags. Therefore the breakdown of this hypotheses signals the limit of applica¬ 
bility of the model. Using this criterion an upper density limit around three 
times the saturation density of symmetric nuclear matter has been found 
[7j. For densities beyond this value the quark-quark interactions through 
the conhnement region should be introduced, so the naive bag picture is not 
sufficient to describe the physical situation in this case. This fact is taken 
into account in the literature using different approaches. For example in [2] 
the quark-quark correlations between bags are introduced in the overlapping 
region. A different approach is given in |H], where the effects of quark-quark 
correlations are assumed to be represented by the exponential dependence of 
the bag constant B on the a meson. 

We propose an alternative viewpoint that preserves the scheme of non¬ 
overlapping bags, and that intends to take care of the strong repulsive com¬ 
ponent of the baryon-baryon interaction that appears as a consequence of 
the internal structure of the particles. This short range repulsion can be 
treated in a simplihed model where baryons are described using a dynamics 
of extended objects. Therefore the fraction of available space is reduced as 
compared to the case of point-like particles. A similar approach has been 
applied to study the phase transition of nuclear matter to the quark-gluon 
plasma imini and in heavy-ion collisions nsnn. 

Since hnite size baryons are assumed as non-overlapping, their motion 
must be restricted to the available space V dehned as cniia 

V' = V (2.19) 

b 


with the total number of baryons of class b inside the volume V, and Vb 
the effective volume per baryon of this class. Hence we conjecture that one 
can renormalize the particle (antiparticle) wave function replacing V for V, 
and thus the effective baryon helds T can be written as 

I a^k, s)u\k, 

k,s 

+ b’’\k,s)v\k,s)e^’^^^>^] (2.20) 

in terms of the Fock space operators a and b, for particle and antiparti¬ 
cle respectively. In this way the hnite size of the baryons is automatically 
accounted for into the held dynamics. 

Eq. ()2.2()j) reinforces the fact that, within this model, the inner and outer 
regions of the bag must be regarded as complementary in order to avoid 
inconsistencies. On the contrary, if these regions can overlap the validity 
of the linear (Eq. (I2.7j) l and non-linear (Eq. (I2.9j) l boundary conditions 
should be revised and additional degrees of freedom, such as exotic multi¬ 
quark states, should be included. Therefore our assumption extends to higher 




densities the quasi-particle picture and it generates further nonlinear baryonic 
couplings. 

For the moment the volumes Vb are parameters associated with the trial 
quantum state of the whole system, and will be determined using a variational 
principle, as it is explained later. 

It is interesting to note that for a mixture of different baryons the excluded 
volume is not exactly the same for all the species ^H]- To simplify the 
discussions, in this paper we neglect these small differences. 

The effective volume per baryon Vb is proportional to the actual baryon 
volume, i.e. for spherical volumes of radius Rb 

47r 

Vb = a—Rb\ (2.21) 

and for sharp rigid spheres a is a real number ranging from 4, in the low 
density limit, to Sa/^/tt, which corresponds to the maximum density allowed 
for non overlapping spheres, in a face centered cubic arrange. Since we wish 
to study the high density regime of homogeneous isotropic matter, we shall 
adopt a = in all our calculations. Thus Vb = 4A/2i?fe^, and the limit 

of validity of the calculations would correspond to densities such that 
the center of mass of baryons are at a distance greater than 2R apart. This 
gives Umax = !{^Rmax^)i where Rmax denotes the biggest radius among all 
present baryonic classes. As we shall see below, with the implementation of 
the procedure outlined in Eq. this limit is never reached in the range 

of densities explored in the present calculations. 

In order to see how excluded volume corrections appear in our approach, 
we shall use the renormalized held of Eq. to calculate the relationship 

among the baryonic densities and the Fermi momenta kb 


C Ic ^ 

/ dx^ < >= (1 - ( 2 . 22 ) 

Jv V 

where ^ —>• V/ (27r^) J dk^ has been used. 

k 

This result is equivalent to Eq. if the factor -d takes the value 

^ = (2.23) 

b 


for hnite baryonic effective volumes Vb ■ In the limit Ub —0 one recovers the 
point-like expressions. 

Eq. fl2.22j) shows that these kind of short-range correlations couple non 
linearly the baryons among themselves, in a density dependent way. Thus 
baryons are considered in this scheme as quasi-particles dressed with these 
corrections. 

The density Ub for a given baryonic species appears on both sides of 
Eq. fl2.22j) . and it is possible to solve it exactly for Ub in terms of all the Fermi 
momenta , namely 


n = 


(1 + E 


ky A 




(2.24) 


b' 







Since depends explicitly upon the baryonic densities, the chemical poten¬ 
tials get an extra term, i.e. 



(2.25) 


(2.26) 


where the energy density en is given by Eq. (l2Tfl) with dehned as in 

En. dT^ . 

The total pressure acquires an additional term AP as compared to the 
pressure of point-like baryons Pq in Eq. (ITTHl) 



(2.27) 


b 


We proceed to give a physical interpretation to the volumes Vf, of the 
baryons, or equivalently their radii Rb- They are considered within this ap¬ 
proach as variational parameters of the trial quantum state of the system. 
Their equilibrium values, at zero temperature, must be determined by min¬ 
imizing the total energy of the whole system, in agreement with nulla. 
This procedure dehnes unambiguously the equilibrium radii of the bags, that 
depend on the baryonic density. 

The new equilibrium condition for the bag consistent with the excluded 
volume corrections introduced in Eq. has the same form as that given 

in Eq. (im . but now ^ = n\/{an^). 

As we shall see in the next sections, as the density grows a state of de- 
conhned quarks becomes more favorable than a system of quarks conhned 
within baryons. Because of this, the concept of excluded volume is meaning¬ 
less when the baryonic phase solves completely. 

In our approach only the baryonic states receive an explicit correction due 
to short range correlations, by the normalization of the helds It must be 
stressed that the meson mean held values are completely determined from the 
baryonic sources, that already include hnite size corrections in d. Leptons 
do not couple to strong interactions, and in this sense they are taken as 
point-like particles. 

To summarize, in the QMC model extended to high densities the hadronic 
matter properties are determined applying the set of equations (ITTIl to (IZTTIl . 
together with (IT^ - dTTTIl for a hxed total baryonic density n, using the value 
of ■§ given in (I2.2d|l . 


3 Quark matter phase transition and the struc¬ 
ture of neutron stars 


Neutron star matter is electrically neutral and it has reached equilibrium 
against /d-decay. The relative abundance of the different baryonic species 












are determined by these conditions. The structure of neutron stars with 
hyperon contributions has been widely studied, for recent investigations see 
for example imni and references therein. In particular in [S] a extension of 
the QMC model is used, that takes into account the density dependence of 
the bag constant B and includes strange mesons. 

In the present work we consider the nucleon duplet (n,p), the A-hyperon, 
the S triplet, the S duplet, and two lepton species, electron and muon. 

In hadronic matter at chemical equilibrium the following relationships are 
fulhlled for the baryonic and leptonic chemical potentials 

fJ’ = /i” + Qfe/i® (if baryon class b is present), (3.1) 

(if muons are present), (3.2) 

with Qh the electric charge in units of the positron charge for the class b of 
baryons, /i^ given by Eq. and = k^ih) = \/mr^ + . 

Mean held equations are solved for hxed total baryonic density n, and 
zero total electric charge density 

n = (3.3) 

b 

0 = (3.4) 

b I 

where the sums run over the baryonic octet and over the two lepton species, 
respectively. 

Once Eqs. (I2.7ll - (l2.9jl have been solved, from Eqs. (I2.1()|1 - (I2.12|) we de¬ 
termine the meson helds (uo, ujq, Bq) and together with Eqs. (I2.15|1 . (I3.1I1 - (I3.4|1 
we get the baryonic and leptonic densities (n^,n^), all the Fermi momenta 
{kb,ki) and chemical potentials (/i^,/i^). 

As the baryonic density increases a phase transition from hadronic to 
quark matter, made up with deconhned quarks, can be reached. Therefore 
the previous set of equations must be modihed to satisfy the new conditions 
of /9-equilibrium in the transition region. Since the baryonic number and 
the electric charge neutrality are always conserved, a smooth crossover be¬ 
tween hadronic and quark matter must be expected because there are two 
conserved charges. In fact, in this coexistence region of mixed hadron and 
quark phases the total baryonic density and the total (zero) electric charge 
are shared between this two phases. 

To determine the composition of the transition region we suppose that 
quarks can exist either in the conhned phase (inside baryons) or as decon¬ 
hned particles. As a hrst approximation we suppose that the deconhned 
phase contains free quarks and leptons, and non-perturbative gluon ehects 
are represented by the bag constant B. 

Under equilibrium conditions, the net hux of quarks through the bag 
surface is zero, and hence the two phases coexist with the only constraint of 
globally conserved electric and baryonic charges. 






We adopt the treatment of (20], as it is appropriate to describe phase 
transitions with more than one conserved charge. Thus, in the coexistence 
phase, the conservation equations (EH and dH are generalized to [20] 


n = {I- x)nH + X^Q, 


0 = + 

b q I 


(3.5) 


where nn = ~ are respectively the hadron and quark 

contribution to the baryon number density; n'^ = Nck ^/is the number 
density of quarks for = 3 colors, kg is the Fermi momentum , and 
is the electric charge in units of the positron charge for quarks of flavor q. 
The quantity y is the volume fraction corresponding to the quark matter 
phase (0 < X < 1). The /^-equilibrium condition for quarks in the deconhned 
phase reads: Hd = + l^ei with /ig = ^rrig + k"^. In the mixed phase 

these relations must be supplemented with //„ = ^fid — and the mechanical 
equilibrium condition Ph = Pq, where Ph,q are the pressures in each phase. 
It must be noted that the case x = 0(1) in Eq. (jd.5|l corresponds to the pure 
hadronic (quark) matter phase instance, the mixed phase is in between. 

The energy density in the mixed phase can be similarly expressed as 
e = (1 — x)^H + X^Q- The quantities eg and Pq are, respectively, the energy 
density and pressure for the deconhned phase which contains free quarks and 
leptons [2D] 


AT r^Q 

B -\ - h / dkk^ Wniq'^ + k"^ 

^ g io 


+ \ [ dkk"^ \/mi^ + k"^, 

(3.6) 

^/Xgu'? + ^/i^n' - eg. 

(3.7) 


q I 


Non-perturbative effects in eg arising from the gluons are represented by 
the bag constant B m- 

Of course, when the pure quark matter phase is reached (x = 1), the 
conditions of chemical equilibrium apply among deconhned quarks, since 
this means that baryons have solved completely. 

The EOS emerging from this calculation can be used to evaluate the prop¬ 
erties of neutron stars. The stellar radius R and mass M are obtained by 
solving the Tolman-Oppenheimer-Volkoh relativistic equations for a spheri¬ 
cally symmetric (non-rotating) neutron star 

, , 2 n K?") + T’(r)] [m{r) + 47rr^P(r)/c^] 

r‘^[l — 2{G / c^)m{r) / r] 

dvrr'^ [€-{t')/(?] dr'. 



dP 

dr 

m{r) 


(3.8) 






Starting from a given valne ec for the central energy density, these eqnations 
are integrated ontward nntil a radius R{ec) is reached for which the pressure 
P is zero, and M = m{R) is dehned. 


To determine more accurately the radius of the star, it becomes necessary 
to use an appropriate EOS for low densities. We have selected the EOS given 
in reference for baryonic densities below O.lno- It is worthful to mention 
that another possibility is to choose the EOS given in which proves to 
be very similar to the former EH- 

The moment of inertia I for a slowly rotating star can be obtained using 

EH EH! 


I 


— [eir)+P{r)]/c^ 

3 Jo — 2{G/c^)m{r)/r 


u}{r) 

~vr 


dr , 


(3.9) 


where O is the uniform angular velocity of the star as seen by a distant inertial 
observer, O -C (c/i?) \f{GjW)MjR. The radial function z/(r) is solution of 
the differential equation 


^ = 2 rr/rO +47rr=^P(r)/c^] 

dr r2[l — 2(G/c^)m(r)/r] ’ 

with the boundary condition v{R) = — ln[l — 2{G/c^)M/R\. The relative 
angular velocity uj{r) measured with respect to the local dragged inertial 
frame, is the solution of 



d 

dr 


4 ■/ 

r j( 


duj{r) 

dr 


+ 4r^a;(r) 


dr 


0 , 


(3.11) 


withj(r) = yi - 2(G/c^) m(r)/r exp [—z/(r)/2] inside the star. Since j(r) = 1 
for r > R, uj{r) has the form a)(r) = r2[l — 2[GI(?') J/r^] outside. 

A useful quantity in astronomical observations is the surface redshift z 
given by 

z = (3.12) 

Numerical values obtained in our approach can be found in Table |21 


4 Numerical Results 

Within the present model the quark masses take the current values = 
rrid = bMeV and rus = ISOMeE. The parameter B is the same for all the 
baryon bags, and we avoid any speculation about its density dependence. 
For a given value of B the set of parameters Zob are adjusted to obtain the 
experimental baryon masses at zero density. 

Since mesons interact directly with quarks, the corresponding meson- 
baryon couplings are related to the quark-meson couplings g'^ (0 = a, a;, p, 
q = u, d) in a. simple way [H] . Denoting as the coupling of the 0-meson to 
6-baryon, 











where is the non-strange quark number inside the baryon b. Thus, for 
given g^, and 5 ^“, the full set of baryon-meson couplings can be determined. 
Their numerical values are obtained by reproducing the symmetric nuclear 
matter properties at saturation, i.e. baryonic density, binding energy and 
symmetry energy 


no = 0.15/m 

Eb = (e/'n.)o — Mc^ = —16MeV, 

where t = [rin — np)/n and M = 938.92MeV/c^ is the averaged free nucleon 
rest mass. 

As expected, the values of the couplings are sensitive to either the inclu¬ 
sion or not of the excluded volume corrections. Both instances are considered 
in Table lU 

For practical applications, we have chosen the values = 169.93,187.83 
and 210.85Meld that will be denoted respectively as (a), (b) and (c) in the 
following. They yield a proton bag radius Rp = 0.8, 0.7 and 0.6 fm respec¬ 
tively (see Table [T}. The cases with the excluded volume corrections (CC) 
have been compared with calculations without them (NC). We remark that 
in the cases labeled NC we take ^ = 'd = 1, and = 0 for all baryons. 

The same B is used consistently in the QMC and deconhned quark de¬ 
scriptions. 

In all cases we start in the hadronic phase increasing the baryonic density 
until the coexistence conditions with quark matter are satished. The excep¬ 
tion here is the case (c)-NC that stays always in the hadronic phase, at least 
within the range of validity of the non-overlapping bags assumption. 

This phase transition occurs at critical densities rid which are shown in Ta¬ 
ble n it can be seen that Ud increases with B. The lowest value Ud — Uq 
corresponds to B^/'^ = 169.93MeId, and therefore we interpret that it is not 
suitable for a true physical consideration and has been kept only for illustra¬ 
tive purposes. 

For the bag parameters B considered here, the non-overlapping bag hy¬ 
potheses breaks down in the NC approach at limiting densities Umaxi before 
the pure quark plasma state has taken place. These upper limits are dis¬ 
played in the last column of Table ^ Thus, stars with a pure quark core 
could not be described within this framework. 

On the other hand, in the CC instance the upper density threshold ricu for 
the mixed phase is enlarged as B increases. No bag overlapping is produced 
in the calculations using the CC option for densities n < ricu, and for n > ricu 
excluded volume effects are meaningless because only quark matter without 
hadrons can be found. 

In Fig. [U we plot the contribution of the hadronic phase {uh) and the 
free quark phase (ng) to the total baryonic density in the mixed phase, see 
upper line of Eq. (EH). For the sake of completeness we draw a dotted line 
corresponding to matter in a pure (either hadronic or quark) phase, all curves 
above (below) this line represent the contributions nq {uh)- The case (c)-NC 



stays in the pure hadronic phase for all the range of densities studied, and 
therefore it coincides with the dotted line in this plot. It can be seen that 
hh lies always below 5.6 no, which corresponds to the upmost case ((c)-CC)) 
considered here. 

In Fig. 121 we show the meson mean held values in the hadronic and mixed 
phases. Differences between the CC and NC treatments become appreciable 
for n > Uq. As can be seen, the p meson amplitude has a sudden change of 
slope at the phase transition due to the change in the isospin composition of 
the hadronic sector in the neighborhood of rich In the same hgure the factor 
d containing the excluded volume correction is displayed as a function of the 
baryonic density. This factor seems to become almost constant at sufficiently 
high densities. 

The proton bag radius Rp and its effective mass M* can be examined in Fig. 
El A faster radius decrease is observed in the CC instance, as compared with 
the respective NC calculation. This fact explains the absence of overlap and 
the consequent lengthening of the validity range of the CC approach. A drop¬ 
ping of about 15% — 20% in Rp is predicted before the hadronic matter has 
completely disappeared; this strong compression of the baryon bags helps the 
deconhnement mechanism in the CC case. This can be understood because 
at a given density the pressure and the chemical potential of baryons are 
slightly greater when the excluded volume correlations are included, than in 
the NC results. Bearing in mind that at the transition point there is a cross¬ 
ing of the curves representing the pressure in terms of the chemical potential 
for the hadronic and quark phases, this crossing occurs at lower densities for 
the CC than for the NC case. On the other side, due to the slowly decrease 
of the bag radius in the NC approach, the overlapping of bags can happen 
before the phase transition to pure quark matter takes place. 

The effective mass M* exhibits a monotonous decrease as a function of 
the baryonic density, the rate of variation at low densities is attenuated by 
both decreasing B and/or including volume corrections. 

The composition of star matter is depicted in Fig. |2]for cases (b) and (c). 
The onset of the quark phase in the NC instance suppresses the hyperons 
that could be present if this transition would not happened. Otherwise in 
the CC approach new hyperon species appear even in the mixed phase. It 
must be noted the pronounced decrease of the lepton abundance in the mixed 
phase. 

The equation of state for (b) and (c) is represented in Fig. El The energy 
density e in the mixed phase varies in the range 0.3 — l.lGeVfm~^ for (b)- 
CC, whereas it ranges between 0.6 — 2.GeVfm~^ for the set (c)-CC. In the 
CC treatment the pressure shows sudden changes of slope at the extreme 
points of the mixed phase, which are absent in the NC case. 

The results for neutron star masses and radii based on these equations 
of state are displayed in Fig. El They are given in units of the solar mass 
Mq = 1.9889 lO^fcc/. _ 

The masses M, radii R and moments of inertia I for the maximum mass 
star are listed in Table 121 We have found that 1.51 < M/Mq < 1.88, a 
result which is above the experimental lower limit M/Mq = 1.44 accepted 
for binary pulsars. 


For these particular cases we have obtained that the radii R are rather 
insensitive to the star internal structure, ranging between 12 km < -R < 
13 km, being slightly bigger when excluded volume corrections are included. 
Fig. El also shows that the maxima are reached in a plateau, and they are 
enhanced in the CC cases studied here with respect to the NC ones (Table 
E}. It can also be appreciated what happens when excluded volume effects 
are not considered. In fact, the NC EOS (dashed lines) predicts that for 
stars approaching the maximum stable mass M the baryonic bags at the 
center are close to overlap. This is indicated by vertical double bars on 
the corresponding curves. In particular for the case (b)-NC the bags are 
overlapping at the center of M, see Tabled (rimaa:) and Tabled (ricmaa:)- On 
the other hand, when CC corrections are included no overlap is present and 
this is one of the points in favor of the inclusion of the excluded volume 
correlations. 

At this point we can compare our predictions with previous works, as 
for example jH] where the QMC model is used to study hyperonic matter 
and neutron star properties. There the bag constant depends on the scalar 
mesons fields, and the increase of the radius with density leads to an early 
overlap of the bags. This feature contrasts with the monotonous decrease 
of the bag radii we have found. We can also examine the EOS, although 
the transition to quark matter is not considered in [Hj. The model named 
QMCI in this reference and the set (b)-NC predict a similar EOS, and for 
the option (b)-CC we find a slightly stiffer behavior in the pure hadronic 
matter region due to excluded volume corrections. At higher densities the 
appearance of the mixed phase softens significantly the EOS, as it is shown 
in Eig. El With respect to the maximum neutron star mass and size in the 
case (b)-NC (Table Ej) compares fairly well with those given by the model 
QMCI 0. 

The structure of the star, in terms of shells of either pure quark, mixed 
quark-hadron, or pure hadron matter, depends on the selected value of the 
bag parameter B. This can be seen in Fig. |71for the CC (b) and (c) cases, 
where we plot the mass m(r) (left panel) and moment of inertia J(r) (right 
panel) enclosed within a given r, for the maximum star masses (Table El) • 
The mass m(r) is given in Eq. 13.81 and /(r) corresponds to the definition 
given in Eq. 13.91 with the replacement R ^ r m the upper integration limit. 
As indicated in this figure, the case (b)-CC predicts that a mixed phase of 
deconfined quarks and hadrons can be found in a central core with a radius of 
about 7.48fcm, which contains 45% of M and contributes 18% to the moment 
of inertia I of the star. Meanwhile in the (c)-CC case the hadron-quark mixed 
phase extends up to 4.31/cm and it encloses only 12% of M; its contribution 
to / is small, around 1.5%. 

For the sake of completeness, we mention that the option (a)-CC, al¬ 
though not physically reliable, predicts for the star with the maximum mass 
a structure which includes a pure quark core of 6.55km. This core is sur¬ 
rounded by a crust of mixed phase located between 6.55km and S.OAkm, 
and the outer shell contains only hadronic matter. The mass of this star is 
1.42M0 with a radius of 9.28km. 


5 Conclusions 


We have studied the high density regime of matter in a schematic model 
of quarks conhned into bags, the so called QMC model. One of the main 
hypotheses of this model is that bags do not overlap, but this situation is 
reached at relatively low densities of around three times the nuclear mat¬ 
ter saturation density in the fixed B treatment. Two different approaches 
have been given to solve this problem by relaxing the requirement of non¬ 
overlapping bags mm However these procedures could give rise to formal 
inconsistencies such as the non-validity of the boundary conditions at the 
bag surface. This point should be clarihed before practical applications at 
extreme densities can be done. 

In the present work we intend to retain all the formal aspects of the 
original QMC, including at the same time the dynamical effects arising at 
very high densities. We do this by modifying the standard QMC treatment 
in two aspects, in first place we consider an alternative condition for the 
equilibrium of a bag immersed in a dense medium. On the second place 
we introduce in-medium averaged short range correlations among baryons 
by using an excluded volume treatment. Thus we obtain a quasi-particle 
picture of baryons dressed by these strong correlations, that remains valid in 
all the range of the present calculations. 

We studied the equation of state of hadronic matter and the transition to a 
deconhned quark phase. Our approach links coherently the nonperturbative 
QCD effects represented by the bag constant B, for both phases of confined 
and deconhned quarks. The excluded volume correction becomes effective 
at intermediate densities, contributing to the onset of a quark matter phase. 
The non-overlapping bags hypotheses is verihed over all the range of densities 
for which hadrons are the relevant degrees of freedom. 

We have found that for star matter a mixed hadron-quark phase precedes 
the pure quark matter, in agreement with previous studies. It is also found 
sudden changes in the compressibility at the extreme points of this mixed 
phase, as well as abrupt changes in the density rate of growth of the vector 
iso-vector meson and in the lepton concentration. 

The effects of the hnite baryon volume correlations on the neutron star 
structure have been examined, we have found that the maximum star mass is 
enhanced by both increasing B and introducing excluded volume corrections. 
The CC approach leads to a small increase of the star radius at the maximum 
mass, as compared to the respective NC case. 

We have demonstrated that the CC treatment enlarges the range of appli¬ 
cability of the QMC model, which can be extended to reach the pure quark 
matter phase. In effect, in our approach this transition takes place at bary- 
onic densities where excluded volume corrections are signihcative, and they 
should be included to describe properly the structure of massive systems, like 
neutron stars. 

Our results are based on a schematic model, however it seems to include 
the main ingredients for a reliable qualitative description of the high density 
regime of matter. Further rehnements can be introduced to give more elabo¬ 
rate predictions, such an improved treatment of the quark matter phase and 


the consideration of exotic mnltiqnark states, which will be object of fntnre 
stndies. 
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Rp[fm] 

Case 

y<j 

yuj 


Ud/no 

‘^cuj ^0 

’^raaxj "^0 

169.93 

0.8 

NC 

CC 

5.747 

4.678 

2.756 

1.701 

8.668 

7.600 

H 

3.98 

3.40 

187.83 

0.7 

NC 

CC 

5.858 

5.311 

2.872 

2.372 

8.582 

7.948 

2.10 

1.87 

6.58 

6.40 

210.85 

0.6 

NC 

CC 

5.993 

5.702 

3.007 

2.747 

8.523 

8.151 

3.91 

10.17 

6.20 


Table 1: The quark-meson couplings for each of the three values of the bag constant 
B used in our calculations, and for each of the approaches with excluded volume correction 
(CC) and without it (NC). The following two columns show the lower (rici) and upper {ricu) 
densities of the mixed hadron-quark phase, and in the last column the limiting density 
Umax for the NC cases are displayed. 



Case 

JAIMq 

R[km 


max / ^0 

z 

I[lR^^kg krri^] 

187.83 

NC 

1.506 

11.70 

6.886 

0.2701 

137.8 


CC 

1.592 

13.02 

5.148 

0.2513 

184.1 

210.85 

NC 

1.672 

12.38 

5.800 

0.2897 

181.5 


CC 

1.879 

12.69 

5.507 

0.3330 

227.1 


Table 2: Neutron star properties for each of the bag constants used in our calculations 
and for each of the approaches: with excluded volume correction (CC) and without it 
(NC). The star mass M relative to the sun mass, the star radius R, the central baryonic 
density {ricmax), the surface redshift z and the moment of inertia / (for a slowly rotating 
star), all corresponding to the star with maximum mass. 
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Figure 1: Hadronic {uh) and quark (ng) contributions to the total baryonic density in 
the mixed phase. The results with (without) excluded volume correction are represented 
with solid (dashed) lines. The different bag constant values = 169.9,187.8, and 210.8 
MeV are distinguished with the labels a, b and c respectively. For comparison we have 
also drawn a curve corresponding to a pure phase (dotted line). All curves above (below) 
this line represent the contributions uq {tiu). The case (c)-NC stays in the pure hadronic 
phase for all the range of densities studied, and therefore coincides with the dotted line in 
this graph. The CC curves are plotted up to the upper density threshold Ucu ■ Numerical 
values for the density threshold Uci and ricu can be seen in Tabled 
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Figure 2: Meson mean field values a,ui and bo as functions of the baryonic density 
relative to the saturation nuclear density uq. The line and label conventions are the same 
as in Fig. ^ In th® right lower corner the volume correction factor is displayed in terms 
of the baryonic density. The double vertical bars indicate the limit of the non overlap 
assumption for the NC cases. 
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Figure 3: The in-medium proton bag radius Rp and proton mass M* relative to its 
vacuum values. The line and label convention are the same as in Fig. |21 




Figure 4: The baryon and quark composition of the star matter in terms of the baryonic 
density for the bag constant values = 187.8 and 210.85 MeV, for the NC and CC 
cases. The line convention is explained in each panel. 













Figure 5: The pressure P of star matter in terms of the total energy density e = E/V, 
for the bag constant values = 187.8 and 210.85 MeV, for the NC (dashed lines) and 
CC (solid lines) cases. For the last instance the pressure has discontinuous derivatives at 
rici and ricu (see Table ^1. 
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Figure 6: The gravitational star mass M/Mq in terms of the central pressure Pc (left 
panel) and M/Mq in terms of the star radius R (right panel). It can be appreciated that 
excluded volume corrections enhance both M and R, CC cases (solid lines). Numerical 
values can be examined in Table |21 The double vertical bars have the same meaning as in 
Fig. [3 indicating that in the NC cases (dashed lines) the bags at the center of the star 
are close to overlap for masses near the maximum. 







Figure 7: The star mass (moment of inertia) enclosed within the spherical region of 
radius r (see Eqs. (Tnni and (IT9ll l is shown in the left (right) panel for the sets (b) (solid 
line) and (c) (dashed line) within the CC treatment. These results correspond to the 
maximum star mass M, as it is listed in Table |21 The mixed phase threshold is marked 
with a vertical arrow. 




